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Abstract. The linear and nonlinear properties of large amplitude electron- acoustic waves 
are investigated in a magnetized plasma comprising two distinct electron populations (hot 
and cold) and immobile ions. The hot electrons are assumed to be in a non-Maxwellian 
state, characterized by an excess of superthermal particles, here modelled by a kappa-type 
long-tailed distribution function. Waves are assumed to propagate obliquely to the ambient 
magnetic field. Two types of electrostatic modes are shown to exist in the linear regime, and 
their properties are briefly analyzed. A nonlinear pseudopotential type analysis reveals the 
existence of large amplitude electrostatic solitary waves and allows for an investigation of their 
propagation characteristics and existence domain, in terms of the soliton speed (Mach number) . 
The effects of the key plasma configuration parameters, namely, the superthermality index and 
the cold electron density, on the soliton characteristics and existence domain, are studied. The 
role of obliqueness and magnetic field are discussed. 
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1. Introduction 

A high frequency electrostatic mode exists in plasmas containing two distinct temperature 
electrons (here referred to as "hot" and "cold" electrons) and ions [Bfl]. For these electron- 
acoustic waves (EAWs) to be sustained, the restoring force comes from the pressure of the hot 
electron component and the inertia is provided by the cold electrons. The electron-acoustic 
(EA) wave (EAW) frequency is much higher than the ion plasma frequency, therefore the ion 
species may be considered as a neutralizing background, which does not influence the dynamics. 
The phase speed v p h of EAWs is much lower than the hot electron thermal speed Vth,h and much 
higher than the cold electrons thermal speed v t h, c , i-e., v t h, c "C v p h <C v t h,h, in order to avoid 
wave damping due to resonance with either the hot or the cold electrons. Here, the indices h 

1/2 

and c refer to the hot and the cold electron, respectively. Furthermore, since v t h,h/c ~ T h ', c , the 
hot electron temperature should be much greater than T c ; ab initio kinetic theory studies 
have shown that EAW can exist only when Th/T c 3> 10 and when the cold electron component 
corresponds to 20% to 80% of the total electron population [2H6]. 

Many theoretical investigations have been carried out since the original work of 
Watanabe pQ, where the existence of EAWs was first predicted [2H7J. The occurrence of EAWs 
is now confirmed by FAST satellite observations associated with auroral density cavities [8] 
and also in the mid-altitude auroral zone [9] (the interested reader is also referred to the 
references cited in Refs. [8] and [9]). Moreover, laser plasma experiments [10] have traced 
the signature of EAWs. The nonlinear propagation of EAWs has been investigated by many 
researchers, both in unmagnetized |TTHT5| and in magnetized plasmas [ToT - flB"] . The EA solitary 
wave characteristics in unmagnetized plasma were studied in Ref. |19j . where the existence of 
negative potential electrostatic structures associated with cold electron density compression 
regions was established. On the other hand, magnetized plasma have been shown to sustain 
EA solitary waves, in the form of pulse-shaped small- amplitude solitary waves in the long- 
wavelength limit [T7J[T8]. The evolution of modulated wavepacket propagation into envelope 
solitons due to oblique perturbations of the wave amplitude has also been demonstrated [16J in 
a Maxwellian plasma. Magnetized as well as unmagnetized ions with Maxwellian hot electrons 
were considered in Ref. [17] . where the existence of negative polarity solitary waves was shown, 
when a drift was present. 

The properties of large amplitude EA solitary waves in ummagnetized nonthermal plasmas 
have been investigated in Ref. [15], where nonthermality was introduced via a Cairns type 
non- Maxwellian distribution of the hot electrons which, interestingly enough, account for 
negative potential solitary structures. An analogous study has been presented in Ref. [13j, 
where solitary wave occurrence in unmagnetized plasma was considered via a pseudopotential 
method, with kappa-distributed hot electrons. The existence of positive polarity EA solitons 
was predicted [18] in a nonthermal plasma by introducing a vortex type hot electron distribution 
in a magnetized plasma. 

Abundant space plasma observations [20 -22] and experiments [71 1231425] have reported 
the presence of highly energetic particles. The observed velocity distribution can be fitted 
by a n — type (generalized Lorentzian) function [71 126T - I28] . where k defines the strength of 
superthermality; the lower the value of n, the higher the degree of excess superthermality, 
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or in other words, the deviation from the Maxwellian behavior. Refs. [5H7J, on the one hand, 
and [13], on the other, have studied the effects of ^-distributed hot electrons on linear and 
nonlinear electron acoustic waves in unmagnetized plasmas, respectively. 

Our main aim in this article is to investigate the characteristics of large amplitude 
electrostatic solitary waves of electron-acoustic type in magnetized plasmas characterized by 
the presence of an excess population in the superthermal hot electron component. As a 
hypothesis, the hot electrons are considered to be k— distributed. The manuscript at hand 
is arranged in the following manner. The basic formalism, relying on a multi-fluid model, 
is presented in Section [2J The characteristics of the linear waves are analyzed in Sec. El 
A fully nonlinear pseudopotential-type approach is presented in Sec. H] to study arbitrary 
amplitude EA solitary waves. The existence domains for EAW propagation are delineated 
in Sec. [5] and their parametric dependence on the relevant plasma parameters is analyzed. The 
soliton characteristics are investigated in Sec. EJ Our findings are summarized in the concluding 
Section [7J 



2. Basic formalism 



We consider a three component plasma model consisting of two electron components (referred 
to as "cold" and "hot") and ions. The ambient magnetic field Bo(= Bqz) is assumed stationary, 
pointing along the z— axis. The direction of wave propagation, at an angle 9 to B , together 
with the magnetic field direction, defines the x — z plane. We also introduce the direction 
cosines a and 7 = cos 9 with respect to the x— and z — axis, respectively. 

The inertia of the cold electrons is retained via a fluid model, that satisfies v t h, c v p h <C 
Vth,h- The (inertialess) hot electrons have a non-Maxwellian character, described by a k— type 
velocity distribution, that leads to a hot electron number density of the form [28] 

r e $ 1 ' K+1/2 

n h = n h0 1- — 3 , (1) 

where n^o represents the hot electron number density at equilibrium, ks is the Boltzmann 
constant and Th is the characteristic "temperature" of the hot electrons [28]. On the timescale 
of interest, the ions may be considered as stationary. 

The dynamics of the cold electron component are described by the following (normalized) 
equations 



dn 

~dt 
du 



+ V-(nu) = 0, (2) 
+ (u • V)u - V<f) - Q c {u x z) = , (3) 

P(n-l)+(l-^Y K+1 -1, (4) 



K 2 



where the cold electron number density n c , velocity u c , and electrostatic potential $ are 
normalized as n = n c /n c0 , u = u c /v s and = $/3>o, respectively. Here, Q c = eB /(m e w ?)ft ) = 
Bq / y/ A.imhom e is the (electron) cyclotron frequency normalized to the (hot electron) plasma 
frequency co p h = (Aim^e 2 /me) 1 / 2 . Space and time variables are scaled by the hot electron 
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Debye length X Dh = (/csT/ i /47m/ l0 e 2 ) 1//2 , and the inverse hot electron plasma frequency, uT^ 1 , 
respectively. Finally, the electric potential $ is scaled as = $/$o = e^/ksTh- The 
characteristic speed scale is the hot electron thermal speed v s = {kBTh/m e ) l l 2 . We have 
defined the parameter (3 = n c0 /nho, denoting the cold-to-hot electron population ratio. 

3. Linear waves 

In the small amplitude limit, we may consider perturbations varying as e *( k ' r -^*) ; an d thus derive 
a dispersion relation, relating the wave frequency u to the wave number k. By linearizing the 
dimensionless evolution equations fl2]) - (|4]), we obtain the dispersion relation in the form 

AtfCu 2 cos 2 8\ 1/2 



1 



^CA = 2 ( n c + W 



~,2 



1± 1 



~,2\2 



(5) 



where 



_ 2 pk 



with c K = (k — |)/(ac — |). Note that Eq. (jHJ) essentially provides the wave frequency in an 
unmagnetised plasma [12j[T6]; cf. the results in Ref. [13]. The parameter (3 appears due to the 
fact that only a fraction of the electrons are taking part in the shielding effect. 

The plasma model described above therefore supports two linear modes with frequencies 
ujc and ua- The upper/lower sign in (JSJ) corresponds to the index C/A, respectively, on the 
left-hand side of Eq. (jSJ). Upon restoring dimensions, it is straightforward to show that in the 
Maxwellian limit k — > oo, the dispersion relation (jSJ) is in agreement with the dispersion relation 
given in Ref. [29]. The wave with frequency uc (fast mode) represents electron cyclotron-like 
(EC) waves (ECWs), that are modified by acoustic effects, while the wave with frequency oja 
(slow mode) represents obliquely propagating EAWs. Both modes are strongly dependent on 
the superthermal character of the plasma, and so is the charge screening mechanism, as denoted 
by the appearance of the factor c K in the denominator in (jB]). The phase speed of the obliquely 
propagating EA mode [lower sign in Eq. (jS])] in the limit k 2 c K and Hi 2 VL 2 reads: 



v^=cos9^J7 K . (7) 

In the vanishing magnetic field limit, the dispersion relation ^5} trivially recovers 
unmagnetized high-frequency EAWs [given by relation (JSJ)], which agrees with Refs. [12] 
and [16j. Naturally, the same dispersion relation is recovered upon setting 6 — > in (JSJ), 
for parallel propagation (recall that the Lorentz force has no component in the direction of the 
magnetic field). In this case, in the long wavelength limit, i.e., for k 2 c K , the wave frequency 



becomes cue — > Jj3/c K , so that the phase velocity is given by 



v P h = \/P/c K , (8) 

in agreement with previous results in Ref. [12] (considering c K = C\ therein) and also with 
Ref. [16] in the Maxwellian case (recovered upon setting c K — 1). A detailed analysis of this 
mode can be found in Refs. [T2] and [TB] • As we will see later in Section El the phase speed v p h 
provides the velocity threshold for (superacoustic) EA solitary excitations to exist, and thus 
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corresponds to the real sound (acoustic) speed in the plasma, as will be discussed below. It 
should be added, for rigor, that we need to consider a phase speed much lower than the hot 
electron thermal speed, in order to validate the fluid description [lj for cold electrons (i.e., 
Vth,c *C v p h <C v t h,h)- The requirement ensuring the validity of our model thus translates as: 

c K K-3/2 

We therefore consider values of the cold-to-hot electron density ratio f3 obeying Eq. (j5J), for 
the respective values of the superthermality parameter k, both in our linear and nonlinear 
analyses. For example, for k = 3, 5 and 100, the upper limit ([9]) implies that f3 <C 1.67, 1.29 
and 1, respectively. 

It may be instructive to analyze from first principles the dependence of electron-acoustic 
waves on various relevant plasma model parameters. Linear wave characteristics depend on: 

- electron superthermality (quantified via c K ): the lower the k value the stronger the 
superthermality (the Maxwellian case is recovered for c K — > 1, i.e., k — > oo); 

- the cold-to-hot electron concentration via (3; in fact, (3 should be between 0.25 and 4 
approximately, in order for Landau damping to be minimized (3HE], the upper limit being 
further reduced by expression ([9]); 

- the magnetic field, via the normalized hot electron cyclotron frequency Q c ; and 
— obliqueness (via 6). 




Figure 1. (Color online) Effect of propagation angle.The linear dispersion relation (Eq. [5]) is 
depicted for obliquely propagating electrostatic waves in a magnetized plasma, for j3 = 0.5, k = 
3 and fl c — 0.5. The upper (three) curves for ujc and the lower (three) curves for uja correspond 
to the fast (electron cyclotron-like) and slow (electron acoustic) mode, respectively. The curves 
have been calculated for different angles with respect to the external magnetic field Bq: 6 = 20° 
(solid curve); = 25° (dashed curve); 9 = 30° (dot-dashed curve). 

Effect of oblique wave propagation. Fig. [TJ shows the effect of propagation angle on the 
linear properties of the electrostatic modes found above: the upper branch corresponds to 
obliquely propagating electron cyclotron-like waves, while the lower branch is for obliquely 
propagating EAWs. Increasing the angle with respect to the magnetic field leads to a decrease 
in the frequency u of the lower (EA) mode, and a frequency increase in the electron cyclotron 
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mode, along with the analogous modification in the phase speeds for both modes. This is 
related, effectively, to the fact that the phase speed of the slow mode varies as cos 9, while the 
restoring force, and hence the phase speed, for the cyclotron modes varies as sin 9. Interestingly, 
a similar trend, increasing the frequency gap between the two modes, was found for oblique 
ion-acoustic waves in Ref. [30] (see Fig. 3 therein). 



0.8 r 




Figure 2. (Color online) Effect of superthermality. The linear dispersion relation for 
electrostatic waves propagating at an angle 9 = 30° with an external magnetic field, for (3 = 0.5, 
f2 c = 0.5. The upper (three) curves and the lower (three) curves depict wc and uja, respectively. 
Different values of k have been considered: n — 3 (solid curve); k = 5 (dashed curve); and 
k = 100 (dot-dashed curve). 

Effect of superthermality. The effects of superthermality on obliquely propagating EA and 
electron cyclotron waves are depicted in Fig. [2j 

It is seen that the more strongly non-Maxwellian distributions (lower k, increased excess 
superthermal particles) lead to both modes having lower frequency and phase velocity. 

Effect of cold-to-hot electron density ratio. The dependence of the frequency on the cold- 
to-hot electron concentration ratio is depicted in Fig. |3j The frequency for both slow and fast 
modes increases with (3, since the latter increases the inertia and decreases the restoring force 
underpinning the oscillations. 

Magnetic field effect. The quantitative influence of the magnetic field (expressed by Q c ) on 
the two modes is shown in Fig. HI A stronger external magnetic field increases the frequency gap 
between the two modes (slow and fast dispersion curves in the diagrams), as shown graphically 
in Fig. HI For normalized k > 1, a significant increase in the EA phase speed is seen to occur 
for stronger magnetic field. 

For perpendicular wave propagation, one may consider 9 = tt/2 in Eq. ([5]). The dispersion 
relation for the upper (cyclotron) mode thus becomes 

u 2 c = Q 2 c + u 2 , (10) 

or, in the long wavelength approximation, lo^ = Vt 2 c + On the other hand, the lower 

(acoustic) mode vanishes in this case. 

Concluding, the cyclotron-like mode tends to the wave frequency Cj for parallel propagation, 
which is, in fact, the wave frequency for the EA modes in an unmagnetised plasma. On the 
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Figure 3. (Color online) Effect of the cold-to-hot electron density ratio f3. The linear dispersion 
relation for electrostatic waves propagating at 30° to an external magnetic field, with k = 3, 
Cl c = 0.5. The upper (three) curves and the lower (three) curves depict toe and uja, respectively. 
Different values of f3 have been considered: /3 = 0.5 (solid curve); (3 = 1 (dashed curve); and 
= 1.5 (dot-dashed curve). 




0.5 1.0 1.5 2.0 2.5 3.0 



Figure 4. (Color online) Effect of magnetic field strength. The linear dispersion relation is 
depicted for electrostatic waves propagating at 30° with respect to the external magnetic held, 
for k = 3, fi c = 0.5 and j5 — 0.5. The upper (three) curves represent up, while the lower 
(three) curves represent uja- Different values of il c have been considered: 57 c = 0.3 (solid 
curve); fi c = 0.4 (dashed curve); and il c = 0.5 (dot-dashed curve). 



other hand, the dispersion law recovers y Q 2 . + co 2 for perpendicular propagation. The frequency 
gap separating the two modes at k = is equal to Q c , but varies with k for k > 1. 

4. Pseudopotential formalism for oblique electron-acoustic solitary waves 

We are now interested in large amplitude solitary structures. These may be analytically 
studied via the (so called Sagdeev-type) pseudopotential formalism, first introduced in Ref . [31] . 
Anticipating stationary profile solitary waves moving at a constant speed, all state variables 
are assumed to depend on a single variable, £ = ax + 72 — Mt, where M is the normalized 
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solitary wave propagation velocity scaled by a fixed sound speed (often referred to as the "Mach 
number"), and a and 7 are the direction cosines along the x and z directions, respectively (thus, 
a 2 + 7 2 = 1). We note, for rigor, that the obliqueness (measured via the angle 9 formed by 
the ambient magnetic field Bq and the wave propagation vector k) must be small, in order 
to ensure a negligible current, which is a requirement for the electrostatic approximation to 
be valid for large deviations from equilibrium [32]. This is assumed throughout the nonlinear 
analysis below. 

In order to solve for the anticipated solitary wave solutions, we integrate the continuity 
and momentum equation by considering vanishing boundary conditions at £ — > 00. The quasi- 
neutrality hypothesis (plasma approximation) is adopted here, for analytical tractability, i.e., 
we no longer allow for space-charge effects through the use of Poisson's equation fl3j), but 
instead equate the perturbed charge densities through Eq. (ITT]). Accordingly, the normalized 
cold electron density is expressed as 



K § 



'IV. 



A tedious algebraic manipulation, described in Appendix A, leads us to an energy type integral 
for the electrostatic potential 0, in the form: 

2 




M,K,p) = 0. 



(12) 



This is formally analogous to a classical mechanical problem of motion in a potential field, 
Eq. (|12p representing the energy balance equation for a particle of unit mass located in a position 
(J) at time £ and moving at speed d(f)/df;. The potential energy is given by a pseudopotential 
\I/(0; M, k, (3, 7) which can be expressed as 

M 2 (k- 1/2) rl 1 / rh x-K+i/2 n -3 



/3(k - 3/2) 



1 1 , 



K - 3/2 



x ( | , - K - 1/2 _ . r 2 x [t , 7 2 (1 + W 7 2 (1 + /3)V 

K-3/2J J 1/3 M 2 (3 2 2M 2 fi 2 

7 2 (l + /3) k-3/2 ^ 7 2 + M 2 /3]/ 6 \ 



M 2 (3 2 K-5/2 
7 2 (l + /3)/t-3/2 
M 2 /3 2 k-5/2 

M 2 f3 + 7 2 



M 2 /3 2 
K+5/2 



k-3/2 



k-3/2 
1 



+ 



7 



2M 2 (3 2 



k-3/2 



-2k+3 



k-3/2 



2(1 + 0) 



1 



K-3/2 



1/2 



- + 



K-3/2 



M 2 f3 + 7 2 (1 + 0) ( l ~ ^72) " 2 ( X + Z 3 ) ( X " .^72) 



2(1 + /3) 2 



K-3/2 



1/2 



/?) 1 



K+1/2 



K~\ 2 



K-3/2 



4 ( K _ 2 )«]0 + (4k - 2)0 2 + 4(1 + - 3/2) 2 (l - ^72) 



K+1/2 



^72) -(1 + /3) (1-^72) 



K+1/2 



Obliquely propagating solitary waves in magnetized kappa- distributed plasmas 9 

x 7 2 M 2 7 2 7 ^( 2 + /3)-l-/3 1 

4(k - 3/2) 2 2(1 + /3) 2 2M 2 /3 2 /5(1 + /5) J ' 1 ' 

We may now study the nonlinear dynamics of EAWs by investigating Eq. (TT5]) in terms of 
different plasma parameters within our model. 

For the sake of reference, we may add that for a Maxwellian plasma, i.e., considering the 
limit k — > oo, expression (fl3|) reduces to 

n 2 c r e *M 2 /3 2 r 2 r m 2 m 2 /3 2 



* Maas (<j>;M,P,>y) = + , ^ , ^ x T^T^T? + 



2 [ (e*-l-/3) 3 J l(l + /3) 2 (e*-l-0) 
+ M 2 /3(2 + g) - g(l + 2/3) 7 2 + 7V 2 + 7 2 (1 - e<^ + 0) 2 



\2 



(1 + /3) 2 M 2 M 2 /3 



2 



+ 7 2 



I >^ 2 + 3/3 



et-l-0j (1 + /3) 
2(e*- 1 -0)(M 2 + 7 V 



^ 2/3[M 2 - 7 2 (1 



1-/3 

M 2 /3 J • (14) 

We may remark, not without surprise, that if we set Q c = 0, expression (1T31) reduces to nil, viz., 
^(0) = 0. This fact shows that our model fails to describe solitary waves in an unmagnetized 
plasma, a limitation imposed by the fact of having adopted the quasineutrality hypothesis in 
the algebra. 

4-1- Approximate analysis for small amplitude solitary waves 

So as to obtain some insight into the behaviour of the solitary waves, it is useful to consider first 
a small amplitude approximation. In this limit < 1), one may expand the pseudopotential 
\l/(0) in Eq. ( flBl) as a McLaurin series near </> 0. A brief algebraic manipulation then leads to 

^(0;M, K ,/3, 7 )^^ - - 



2 M 2 ( M 2 — Mf 



f 157 2 + /3(47 2 -l)(2^ + l) 



6(Mf-M 2 ) 2 l 2k -1 

37 2 M| 3M 2 /3 - 4 + 2(4 + /3)/c 



(15) 



M 2 M| 2k - 1 

where we have defined the characteristic speed M2 = [/3(2k — 3)/(2« — 1)] 1//2 , which we recognize 
as the phase speed in an unmagnetized plasma, see in Eq. (JED. Eq. ( fl5T) describes weakly 
nonlinear waves at a speed M (albeit within a range, to be determined below). 

Taking the small amplitude analysis an order higher, we consider the weakly superacoustic 
range, traditionally described by a Korteweg-de Vries (KdV) type approach |33j. The true 
acoustic speed, which provides the lower bound of the Mach number range of permitted values, 
may be found by considering the roots, for a pseudopotential ^((f>), of the equation \l/"(0) = 0. 
Using the above pseudopotential (Eq. [IS]), one finds from ^'/(O) = that the acoustic speed 
for oblique (linear) electrostatic waves is given, in our case, by Mi = I7IM2 = M2COS6: see 
Eq. (JED- 



Obliquely propagating solitary waves in magnetized kappa- distributed plasmas 



10 



Upon setting M ~ Mi+v = \ r y\M2+v, where v <C M\, and linearizing in v, expression ( TT5|) 
is cast in the form: 

W^'^'^^-^Kl-T 2 )^ + 6M 2 4 (l- 7 2 ) 2 l 2k -1 

o o/3 -4 + 2(4 + /3)k1 , , . 

+ 15 7 2 -3 - 3 7 2 ^ 16 

2k — 1 J 

The smallest non-zero root that annuls the pseudopotential yields the amplitude of the 
solitary waves, m , i.e., ^(^m! v, k, (3, 7) = provides the amplitude. This leads to 

\<f> m \ttf{K,(3 n )v, (17) 

where 

6M 2 (l- 7 2 ) 



/(«,& 7) 



x 



7 

2 /3(4 7 2 - 1)(2k+1) 4 + 2(4 + /?)^ 

15 7 — 3 H o 7 



2k -1 2k -1 

and 5 ~ M — M\. It is clear from Eq. (fT7|) that |0 m | oc /(k, /3, 7) [or |0 m | oc F(k,(3, 7)], 
where the proportionality constant takes the form M — Mi or [M/Mi — 1 with F(k,/3, 7) = 

|7|VW/(K,A7)]. 

We shall return to this approximate solution for small amplitude solitons later in the paper. 
5. Existence domain for EA solitons 

We shall now investigate the existence conditions and domains for arbitrary amplitude electron- 
acoustic solitary structures, based on the full nonlinear pseudopotential, Eq. ( fTBl . 

First, let us point out that the pseudopotential given by our plasma model in Eq. ( IT3l) 
satisfies indeed the condition ty(<fi = 0) = ^-|^=o = 0, since = defines the equilibrium 
state. Furthermore, we recall that the maximum (absolute) value attained by the electrostatic 
potential 4> m is given by the root of the pseudopotential ^(4>) at m ^ 0. 

The convexity requirement ^^|^ =0 < (also referred to as the soliton existence 
condition [M]) must be imposed, in order for the origin to correspond to a local maximum 
of the pseudopotential function \I>. Expression (TT3"1) thus leads to the analytical requirement 



#% =0 c M 2 (M 2 -/3/c K ) " liyj 



For the sake of reference and comparison, we may adopt the notation of Ref. [35] (see Eq. 
5 therein) by writing the convexity requirement given by Eq. f fT9|) as 

f(M) = V-^W + J^S- > (20) 

Note that Eqs. (fT9i) and (1201) are only valid for 7 7^ 1; in other words, they are not satisfied for 
parallel propagation, 9 = 0. 

The convexity requirement (1T91) imposes that the soliton speed M value should lie between 
a threshold M 1 and an upper bound M 2 , as 



Mi = |7|y / /3/c^ , M 2 = ^~ K . (21) 
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The two limiting values are simply related via a factor 7 = cos 6, related to the propagation 
direction (with respect to the external magnetic field). Interestingly, the upper limit M 2 is 
essentially the phase speed of electron-acoustic waves in an unmagnetized plasma (see discussion 
in section [3]) and does not depend on the direction of propagation. On the other hand, the 
threshold M\ is the real sound speed for obliquely propagating (linear) waves - cf. (|7j) above - 
thus naturally, solitons occur above this limit. 

Summarizing, EA solitons exist in the velocity range 

Mi < M < M 2 , (22) 

■ <*> 

One notes that, for any values of k and (3, the effect of increasing the value of 7 (i.e., of aligning 
the wave vector to the direction of the magnetic field) is to shrink the existence region for 
EA solitary structures. We add, for rigor, that the existence domain vanishes in the limit of 
parallel propagation (see the definitions of M\ and M 2 , which then coincide for 7 = 1), and the 
analytical model breaks down. 




Figure 5. (Color online) Variation of the lower threshold Mi (lower curves) and the upper 
limit A/ 2 (upper curves) with the superthermality parameter k for 30° obliquely propagating 
waves. The curves correspond to (3 — 0.5 for the solid (red) curves; (3 = 0.7 for the dashed 
(blue) curves; and (3 = 1 for the dot-dashed (green) curves. 

The soliton existence region [Mi, M 2 ] is depicted in Fig. [5] for a representative choice of 
the parameter values, for different values of the cold electron concentration /3. We see that M 
(G [Mi, M 2 ]) varies significantly for k < 6, whereas the interval [Mi, MJ\ remains practically 
constant for k > 6. This is due to the fact that the acoustic limit Mi decreases rapidly as 
k — > 3/2. As expected, the plasma behaves as effectively Maxwellian for large values of k 
(practically, above ~ 10). 
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Interestingly, the existence regions predicted for two different values of (3 may not overlap 
at all (observe Fig. E]to see this, for instance), showing that solitary waves propagating at a 
given speed may exist for one value of but not for another. 

We note here that there are a number of analogies between the oblique propagation 
of electron-acoustic and ion-acoustic solitary waves. Formally, the accessible range of Mach 
numbers naturally takes the same mathematical form as is the case for obliquely propagating 
ion-acoustic solitary waves [30]. At the same time, it is clear that the actual velocities would 
be very different, as the normalizing speed here is the hot electron thermal speed, whereas in 
the analogous case it is the much lower ion-acoustic speed. 

We conclude this section by pointing out that the velocity interval (region between M% and 
M 2 ) for the existence of EA solitary waves leads to M 2 — Mi = (1 — cos9)^J(3/c K . The velocity 
interval is therefore expected to be reduced for lower values of /3, obliquity or superthermality 
parameter k. The upper limit M 2 does not depend on the direction of propagation. 

6. Electron-acoustic soliton characteristics 

We are interested in investigating the features of large-amplitude electron-acoustic solitons 
[corresponding to solutions of Eq. ffl2|) . obtained numerically], in fact focusing on their 
dependence on various plasma parameters (superthermality parameter k, cold electron 
concentration via (3, propagation direction via 9). Numerical integration of Eq. ( TT2|) provides 
us with the information about the electrostatic potential and also the corresponding electric 
field excitation E (= — V0). We restrict ourselves to a region where the linear EA wave may be 
sufficiently weakly damped to be observable, i.e., a region of cold-to-hot electron density ratio 
(3 supported by Eq. ((9]) in order to ensure the validity of the fluid model, as well as respecting 
the limits imposed by Landau damping on the linear wave. 

Superthermality effect (via k). We focus on a low range of values for k, to investigate 
the effect of superthermality on EA solitary structures. The soliton speed M is kept constant 
here, within the accessible range of values. Fig. [6] shows the effect of excess superthermality on 
the soliton characteristics, for a fixed speed value (M = 0.53) representative of the accessible 
velocity region (cf. Fig. [5] for f3 = 0.5). One sees in Fig. |6]^a) that both the root and the depth 
of the pseudopotential well increase with decreasing k. This means that both the amplitude 
and the profile sharpness (maximum electric field) of the negative potential EA solitary waves 
(Fig. [6] [b]) increase with increasing superthermality. This is qualitatively reminiscent of the 
results for IA solitons in Ref. [30] (for magnetized plasma) and [36] (for unmagnetized plasma). 

The threshold above which the solitary waves may occur, namely the true acoustic speed 
Mi, has earlier been identified as the phase speed [cf. Eq. ([7]) above] of obliquely propagating EA 
excitations; this depends on both k and /3. One thus expects that the true Mach number MjM\ 
(ratio of the pulse propagation speed to the true acoustic speed for the plasma composition 
under consideration), which differs for different parameter values (i.e., for different k and /3) 
even for a fixed soliton propagation speed M, will be a meaningful physical quantity, whose 
effect on the soliton characteristics should be investigated. Furthermore, it is known [371138] 
that the pseudopotential function satisfies d^/dM < 0, hence d(j) m /dM > J37J[38] (here <f> m 
is the maximum amplitude of the solitary excitations) for barotropic plasma fluids, suggesting 




Figure 6. (Color online) (a) Variation of the pseudopotential ^((f>) with <f> for different values 
of k, and il c = 0.5, /3 = 0.5, M — 0.53, 8 — 30°. From top to bottom: dot-dashed curve for 
k — 4; dashed curve for k = 3.5 and solid curve for n — 3. (b) The corresponding electrostatic 
potential perturbations, obtained numerically. 



that the soliton amplitude is an increasing function of MjM\. In addition, the amplitude of the 
the Korteweg-de Vries (KdV) solution (valid for small potential amplitudes) has a linear 
dependence on the supersonic velocity increment [i.e., oc (M — M\), or oc M\{M j M\ — 1), having 
adopted an appropriate scaling], as shown in Eq. (fT7|) . 

We have numerically obtained and investigated the variation of the soliton amplitude m 
as a function of M — Mi, for different values of the relevant parameters, in particular, the 
superthermality parameter, k, and the cold-to-hot electron density ratio, [3. The dependence 
of the acoustic speed Mi on (3, k and 7 has been taken into account, for a given plasma 
composition. Fig. [7J shows the soliton amplitude (f) m against M — Mi and against M/Mi for 30° 
propagation (the true Mach number lies in the range [1, 1.1547]; see definition in Eq. ( 123]) ). using 
the values (3 = 0.5 and Q c = 0.5. Both graphs show practically linear variation of the amplitude, 
but with the slope being significantly K-dependent, particularly for small k. The dotted lines 
represent best fits, and in Fig. [7(a) yield fitted slopes of 0.95, 1.11 and 1.30, for the three values 
of k = 3, 5, and 100, respectively, the last value clearly describing a quasi-Maxwellian case. As 
the amplitudes are relatively small, one may expect that the KdV approximate solution should 
provide a reasonable description of the figures. Indeed, substitution of the parameter values in 
Eqs. (I17|I18I) leads to slopes of 1.02, 1.19 and 1.39, respectively. Similar results are found when 
fitting to the slopes in Fig. [7J(b). We note that the KdV expressions over-estimate the observed 
slopes by a few percent. This is understandable, because the actual curves slowly flatten out 
and slip below a true linear form as M — Mi, and hence \4> m \, increase towards the limit of 
applicability of the small amplitude approximation. This leads, for instance, to the observation 
that, for (M — Mi) < 0.06, the fitted (dashed) curve lies marginally below the calculated roots. 

At first sight, the results presented in the Figs. [6] and [7J appear to be contradictory. From 
Fig. El one might conclude that lower k implies stronger solitons, for fixed speed M. On the 
other hand, Fig. [7J suggests exactly the opposite trend, in that lower k implies lower amplitude 
\<f>m\ f° r given M/Mi. The explanation lies in the fact that the former figure relates to a fixed 
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(a) (b) 

Figure 7. (Color online) Variation of the pulse amplitude \<p m \ (absolute value) with (a) the 
super-acoustic excess speed, M — Mi, and (b) the true Mach number M/Mi, for different 
superthermality parameter values k, as shown on the figures. The other parameters are fixed 
at j3 = 0.5, fl c = 0.5 and = 30°. Note that the curves have been truncated at the maximum 
accessible values of M/Mi, namely at = M 2 /M 1 = 1.15466, 1.15469, 1.1547 for k = 3, 5, 100, 
respectively. The dashed lines represent the best fits to obtain values for the slopes. 

speed M, which implies a variable excess speed (or the true Mach number in Fig. [7]) relative to 
the the true acoustic speed for the plasma composition under consideration. Thus the changes 
observed as k is varied relate strongly to the effect of k on Mi. On the other hand, in the 
latter figure, that effect is eliminated and one is left with the changes in amplitude that follow 
on the true variation from the sound speed (or the true Mach number) alone. Therefore, while 
Fig. [6] investigates the soliton characteristics in terms of a constant speed scale, Fig. [7] adopts 
a variable, k— dependent yardstick, related to a fundamental physical unit, the true acoustic 
speed. 

In qualitative terms, smaller values of k do indeed lead to more intense solitary waves 
(potential pulses), for given speed M, as Fig. [6] suggests; on the other hand, the solitons are 
then "less super acoustic," in loose formulation, and thus inevitably of smaller amplitude (and 
this, for given value of MjM\). 

Effect of the cold electron concentration (via (3). We next turn to an analogous study of the 
effects of the cold-to-hot-electron density ratio (j3) for a strongly superthermal plasma (k, = 3). 
We assume parameter values k = 3,M = 0.53, Q c = 0.5,8 = 30°, and vary over a narrow 
range. One notes that, as (3 = n 0c /n h is decreased, the relative inertia (the cold electron 
mass density governed through uq c ) is decreased, while the relative restoring force (the hot 
electron pressure governed through n^) is increased. These two effects both tend to increase 
the frequency of the electrostatic oscillations underlying the EA wave and soliton. Hence, for 
the three cases illustrated in Fig. [BJ the true Mach number MjM\ increases with decreasing /3, 
as M/Mx = 1.075, 1.095, 1.117 for = 0.54, 0.52, 0.50, respectively. As a result, one would 
expect both the soliton amplitude and the maximum soliton electric field to increase. In terms 
of the pseudopotential curves, these effects imply that one would expect an increase in both 
the position of the root and the well-depth of the pseudopotential. We immediately recognise 
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Figure 8. (Color online) Variation of the pseudopotcntial ^>(4>) for different values of j3, with 
cj) for fl c = 0.5, k — 3, M — 0.53, 9 — 30°. From top to bottom: the dot-dashed curve for 
P = 0.54; the dashed for (3 = 0.52 and the solid for /? = 0.50. 
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Figure 9. (Color online) Variation of the pulse amplitude \4> m \ (absolute value) with (a) the 
true Mach number MjM\\ and (b) the super-acoustic excess in speed, M— Mi, for different cold- 
to-hot electron population ratio /3. The other parameters are fixed at k — 3, Q c — 0.5 and 6 = 
30°. Note that we have truncated the curves at M/Mi = M 2 /M\ = 1.15466, 1.15354, 1.15458 
for /3 = 0.5, 1, 1.5, respectively. The dashed lines represent best linear fits to the data. 



In Fig. [9] we next present figures that are analogous to those of Fig. [3, but varying (3 
rather than k. Other parameter values are fixed at k = 3, Q c = 0.5,6 = 30°. It will be seen 
that the values for (5 that are used in Fig. represent cases in which the cold and hot electron 
densities are equal, as well as cases in which the cold and the hot electron densities, respectively, 
dominate: the values used are (3 = 0.5, 1.0, and 1.5. As we have seen above, the sound speed 
(soliton velocity threshold), Mi = [f3(2n — 3) /(2k — l)] 1 ^ 2 - One would thus expect behaviour 
with respect to increasing /3 to be similar to that seen for increasing k in Fig. (EJ). 

The variation of <p m with M/Mi and with M — Mi is explored in Fig. |9j The figure shows 
that with a decreasing cold-to-hot electron density ratio, the true Mach number, and hence 
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the maximum amplitude of the soliton increases (see Fig. [9]). The best linear fit in Fig. E^b) 
yields values for the /3-dependent slope of 0.95, 1.15, and 1.23, respectively, as (3 is increased. 
These values are well reflected by the KdV solution, based on Eqs. (I17|ll8p . which yields slopes 
of 1.02, 1.2 and 1.28, respectively. Again, one finds that the KdV solution over-estimates the 
fitted slopes by a few percent, as discussed above for the case in which k was varied. 

Obliqueness effect (via 9 or j). We have analyzed the effect of obliqueness of propagation 
angle on electrostatic solitary waves, for given values of the other plasma parameters. The 
Sagdeev potential, as well as the electrostatic potential (pulse), are amplified with increasing 9, 
as shown in Fig. [LTJJ In particular, both the root and the depth of the pseudopotential (\l/) well 
are increased, thus leading to larger values of both the maximum derivative d<p/dl; (maximum 
steepness of the soliton profile) and the soliton amplitude \4> m \, as shown in the plot. This 
is intuitively expected. We recall that the critical Mach number threshold Mi decreases with 
increasing obliqueness as Mi oc cos 9. Therefore, the difference M — Mi for fixed M increases 
with 9, resulting in both a deeper Sagdeev pseudopotential and to a larger root, and thus to 
stronger electrostatic excitations. 
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Figure 10. (Color online) Obliqueness effect, via 9 (the angle between wave propagation 
direction and magnetic field) : (a) variation of ^(<j>) with </> for varying 6, and il c — 0.5, 
K = 3, M = 0.5, j3 = 0.5. From top to bottom: solid curve for 30°; dashed curve for 32° 
and dot-dashed curve for 34°. (b) The corresponding solutions for the electrostatic potential 
perturbation, obtained numerically. 

Effect of Magnetic field (via Q c ). The existence domain [Mi, M 2 ] [discussed in Section El 
see Eq. (T2~TT) ] for the existence of EA solitary excitations does not depend on the magnetic field 
Bo (since it does not depend on the value of fl c ). However, the depth of the pseudopotential 
well increases for larger Q c , i.e., for stronger magnetic field B , as shown in Fig. [TTJ The 
corresponding soliton solutions thus yield constant amplitude, but reduced soliton profile width 
(i.e., are more sharply peaked), for larger Q c . 
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Figure 11. (Color online) Magnetic field effect (dependence on J7 C ): (a) variation of &((f>) 
with (j) for k = 3, M = 0.53, f3 = 0.5, 9 = 30°. From top to bottom: fl c — 0.5 (solid curve); 
fl c = 0.6 (dashed curve) and Sl c = 0.7 (dot-dashed curve), (b) The corresponding electrostatic 
pulse excitation, obtained numerically. 



7. Conclusions 

The characteristics of large amplitude obliquely propagating EA solitary waves have been 
investigated in a three component magnetized superthermal plasma. 

We first briefly considered the linear modes and found two distinct types, namely, 
obliquely propagating magnetized electron-acoustic waves and electron cyclotron-like waves. 
The frequency gap between the two modes decreases with obliqueness, i.e., with the angle 
between the wave propagation direction and the direction of the magnetic field. On the 
other hand, a simultaneous decrease of the frequency of both modes is observed for stronger 
superthermality (lower k), i.e., for larger deviation from the Maxwellian, as well as for an 
increase of the cold-to-hot electron population ratio, /3. 

The main emphasis of the paper is the study of nonlinear electrostatic solitary waves, 
using a pseudopotential method, adapted to treat obliqueness effects. For simplicity, we have 
adopted the quasineutrality hypothesis for our analytical study in the nonlinear part. We have 
determined and analyzed the velocity domain where solitary waves may occur, in the form 
of a threshold Mi and an upper bound M 2 , both of which depend on k and on 0. Here, 
Mi = M 2 cos 9, where 9 is the angle is the angle between the direction of propagation and 
the background magnetic field. The relevant plasma configurational parameters (e.g., the 
plasma superthermality, the cold electron concentration, the propagation angle) are seen to 
significantly modify the allowed regions for the linear EAWs to propagate, along with the 
associated characteristics of solitary waves. 

We have shown that the true acoustic speed of electron-acoustic waves in a magnetized 
plasma decreases with stronger plasma superthermality, decreasing cold electron concentration 
and increasing obliqueness of propagation. This variation of the true acoustic speed with 
different plasma parameters therefore modifies the characteristic properties of solitary waves 
(e.g., amplitude and width). 
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For the sake of rigor, we have tested the validity of our oblique model in various special 
cases. Interestingly, the two limiting Mach-number values (Mi,M 2 ) merge in the parallel 
propagation limit (9 — > 0), indicating that the parallel wave propagation in magnetized plasma 
can not be described by this model. Similarly, the model breaks down in the vanishing-magnetic- 
field limit, which may be attributed to the plasma approximation, here adopted for analytical 
tractability. 

Beyond the strict interest in the context of plasma waves, our study is of methodological 
interest, as it adopts an analytical approach not used before (with the sole exception of Ref. [30J, 
where this first appeared). In particular, the fluid-dynamical equations, which generally cannot 
be disentangled into a single pseudo-energy balance equation if oblique wave propagation in 
magnetized plasma is considered, were managed into a convenient form by assuming charge 
neutrality. There are a number of analytical issues thus raised. 

First of all, a comment is in order, regarding our choice of adopting the plasma 
approximation (neutrality hypothesis, NH). As a matter of fact, the question of whether or 
not the NH is appropriate for the modeling of electrostatic excitations appears to be as old as 
plasma modeling itself, and is discussed in early plasma textbooks [30HT] . yet on a qualitative, 
first-principles basis. As a rule of thumb, it was suggested that high-frequency waves should 
not be modeled by making use of the NH. Later literature has attempted to shed some light on 
the problem, with rigorous emphasis on localized nonlinear excitations (solitary waves), where 
the legitimate quantity to deal with is the phase speed, rather than the frequency (notice the 
discussion in Ref. [42J). Furthermore, in the "conventional" version of the pseudopotential 
method, for magnetic- field-aligned solitary waves, one encounters a second upper limit, related 
to the infinite compression point (or limit): this point cannot be accessed via the methodology 
adopted in this paper. This may be attributed to the physical and also analytical constraints 
implied by adopting the plasma approximation (neutrality hypothesis). 

Admittedly, our choice of adopting the NH was at first dictated by algebraic tractability (as 
the 2D fluid equations could not be disentangled into a simple evolution equation by retaining 
the full Poisson dynamics). The algebra leading from (E])-© via the NH to (|T2l)-(fl3l) is presented 
for the first time, to our knowledge. From a methodological point of view, therefore, we find that 
this study exhibits novelty and interest which should attract future attention and discussion, in 
other nonlinear plasma wave related problems. Focusing on the particular area of electrostatic 
solitary waves, our results on the role of non-Maxwellian behaviour in plasmas include a number 
of clear predictions on soliton existence and stability, which may be tested, e.g., against space 
observations, and hopefully confirmed. 
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Appendix: Derivation of the Sagdeev-type pseudo-energy-balance equations 



Defining the moving coordinate £ = ax + 72; — Mt (recall that 7 = cos 6) one can write the 
normalized fluid equations (j2J and ([3]) in the form 



■^dn d{nu x ) d{nu z 



[— M + au x + 7« 2 



-M + au x + 7« 2 



du x dcj) 

1 a — 

d£ d£ 



0. 

I ^^igUy 



dUy 



f2 c Ma; = 



(24) 
(25) 
(26) 
(27) 



(-M + au x + l u z) ~^f ~ 7^7 = 

Anticipating stationary localized excitations (vanishing at infinity), we integrate Eqs. ()2 
and ( 1271) using the boundary conditions n — >■ 1, — > 0, and <fi — > at £ — )■ ±00, to obtain 



Cfttj; + 7^2 



li- 



lt , 



_ _7 /I 

M \(3 

Mi 1 



+ / n 



a 



7 2 /I 



V nj Ma\j3 



— + nd(f) 



(28) 
(29) 

(30) 



where we have taken quasineutrality, i.e., Eq. ( 1TT1) . into account. Substituting Eq. (128"]) into 
Eqs. ( 125|) and ( )26|) . we can write 

— -a— + i2 c ^ = 0, (31) 

n a£ a£ 







Substituting for m x from Eq. ( 130|) . Eq. ( 1521) can be written as 
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a a M 2 aA/3 
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Differentiating Eq. ( 15Tj) with respect to £ and substituting the value of and du y /d^, we get 

F(0) , (34) 



3M 2 /dn\ 2 M 2 d 2 n d 2 <p 
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(35) 



We may now make use of Eq. (lllj) in order to eliminate (the normalised density) n in 
Eq. (154"]) . leading to an expression in terms of (the electrostatic potential) (f> only, 

( K _ 1)2 r l 1 / ^ x _«+!-, -4/ ^ \ -2k- 1 /rlA\ 2 
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Differentiating twice with respect to £, the above equation reduces to 
M 2 d 2 rl 1 A. \- K +5i- 2 dV 
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d 2 ( M 2 r 1 1 , 



d 2 5 
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Differentiating S with respect to £ and squaring, we write 

r -\ 2 / dsh\ 2 
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where 
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Multiplying both sides of Eq. (I3"8"j) by dS/d£, we can write 



Integrating this equation with boundary conditions — » and dip/d^ — )■ 
one can obtain an energy balance equation in the form 

2 

+ *(0; M, «,/3) = 0, 



at £ 




(36) 
(37) 
(38) 
1 . (39) 

(40) 
—> ±oo, 

(41) 



where \I/(0; M, k, /3) defines the Sagdeev-type pseudopotential of the system. 
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